function d = deviation( s1, s2, time )
    pts = 7:7:time;
    d1 = deval(s1,pts);
    d2 = deval(s2,pts);
    
    m = 0;
    
    for i = 1:length(pts)
        d1cur = d1(:,i);
        d2cur = d2(:,i);
        d1total = RNA2RNATotal(d1cur(1), d1cur(10), d1cur(9));
        d2total = RNA2RNATotal(d2cur(1), d2cur(10), d2cur(9));
        m = max( abs( (d1total)-(d2total) ), m );
    end
    
    d = m;
end
